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sensitive maps with angular resolution often exceeding 50 /ias. Using the GMVA, a large sample 
of prominent 7 -ray blazars have been observed approximately 6 monthly from later 2008 until 
now. Combining 3 mm maps from the GMVA with near-in-time 7 mm maps from the VLBA- 
BU-BLAZAR program and 2 cm maps from the MOJAVE program, we determine the sub-pc 
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can be estimated at different locations along the jet under the assumption of equipartition between 
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consistent with the predictions of magnetic jet launching (e.g. via magnetically arrested disks 
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1. Introduction 

Strong magnetic fields have long been considered important in the production of AGN jets 
around super-massive black holes (SMBH) [4], and in the accretion disk [2]. Recently, there has 
been increasing observational evidence for particularly strong magnetic fields - near fheir fheoref- 
ical maximum - under the model of magnetically arrested disks (MAD) [9, 25]. MAD describes 
a scenario where a very strong magnetic field is produced by accreting maffer being dragged in a 
poloidal magnefic field, causing gas fo fall onto the black hole with a higher-than-free-fall velocity 
[3, 18]. This generates a poloidal magnetic field which is fhen fwisfed around fhe black hole, cre¬ 
ating a toroidal magnefic field in the jel [16, 24]. 

The Global mm-VLBI Array (GMVA) operafes af 3 mm and 7 mm and can provide angular 
resolution exceeding 50/fas af 86GHz [11]. By combining GMVA observations af 3 mm with 
near-in-time observations at 7 mm as part of the VLBA-BU-BLAZAR program [12], we are able 
to spectrally decompose the VLBI structure of y-ray bright blazars at different positions along the 
jet. This allows estimates of the magnetic field sfrengfh fo be delermined af Ihese posilions under 
fhe assumplion of equiparlifion of magnetic field and relafivisfic parficle energies. By making 
some furlher assumptions abouf fhe jef and magnefic field geomefry, we can esfimafe fhe disfance 
befween fhe VLBI “core” (fhe mosf upsfream visible pari of fhe jef in VLBI maps) and fhe cenlral 
SMBH. We can fhen exlrapolale fhe magnefic field in fhe parsec-scale jef back fo fhe SMBH. We 
find fhal fhe mm-wave “core” is localed ~ 1 — 3 pc downslream of the jef base and that the magnetic 
field sfrengfh is of fhe order Rapex ~ 5000 — 20000 G. Such magnefic field strengfhs are high bul 
nol inconsisfenf with those predicted by magnetic jet launching scenarios such as MAD. 


2. Observations 

Six blazars are included in this analysis: 0716-1-714, 0836-1-710, 3C 273, BL Lac, 3C 454.3 and 
OJ287. Three millimetre VLBI data were observed during the period 2008.78-2012.35 using the 
Global mm-VLBI Array (GMVA) and CLEAN maps were produced using the Caltech DIFMAP 
package [22]. Data were fringe-fitted and calibrated using standard procedures in AIPS for high 
frequency VLBI data reduction [e.g. 12] with extended procedures written in ParselTongue as de¬ 
scribed by Marti-Vidal et al. [15]. Within AIPS, amplitudes were corrected for system tempera¬ 
tures, sky opacity, gain-elevation curves and then averaged over all IPs to increase SNR. Manual 
phase calibration was performed on the brightest sources and scans within the experiment. Data 
were correlated at the Max-Planck-Institut fiir Radioastronomie in Bonn, Germany. Compliment¬ 
ing GMVA observations are near-in-time observations from the use of 43 GHz VLBA data from the 
VLBA-BU-BLAZAR Monitoring Program VLBA-BU, funded by NASA through the Fermi Guest 
Investigator Program. The VFBA is an instrument of the National Radio Astronomy Observatory. 
The National Radio Astronomy Observatory is a facility of the National Science Foundation oper¬ 
ated by Associated Universities, Inc. The data were reduced using similar methods described by 
Jorstad et al. [12]. 
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3. Analysis 

3.1 “Core” Identification 

A critical part of VLBI analysis is “core” identification, as it is required in order to perform 
kinematic analysis. The VLBI “core” is assumed to be the most upstream visible component in 
a VLBI image and is assumed to be stationary. The “core” is typically identified on fhe basis 
of i) morphology, ii) small size, higher flux densifies and correspondingly higher brighfness fem- 
perafures fhan downsfream componenfs iii) an opfically fhick (inverfed) or flal specfrum and iv) 
increased levels of flux density variafions. Af wavelengfhs longer fhan af millimefre, fhe “core” 
is oflen relafively easy fo idenlify, as if is usually fhe brighfesf componenf by far. However, af 
high frequencies. Ibis can somefimes be difficull, as fhe sub-parsec scale sfrucfure can include 
quasi-sfafionary fealures fhaf are somefimes brighfer and more compacf fhan fhe “core” [12]. Nev- 
erfheless, based on fhe criferia lisfed above, if has been always possible fo confidenlly idenlify fhe 
“core” in our VLBI maps. 

3.2 Magnetic Fields 

As VLBI at 3 mm allows observations at or above the turnover frequency for synchrotron 
radiation, it allows for a novel approach for deriving estimates of the magnetic field in blazars. 
Although at least three spectral points are required to derive a spectrum (and determine with confi¬ 
dence the turnover frequency), the combined use of 3 and 7 mm VLBI maps allows for limits to be 
computed for the magnetic field as a function of distance down the jet. In the future, the addition 
of 1 mm and 2 cm MOJAVE data will allow us to determine the turnover frequency more robustly. 
For this analysis, only components that could be fitted to a single, non-delta component at both 
frequencies were used. 

Jet emission is said to be in equipartition when relativistic particle and magnetic energies are 
equal. If we assume that the jet is in equipartition, we can compute an estimate of the minimum 
magnetic field strength can hence be calculated from (e.g. [1]): 

Bequi = 5.37 X 10l2(5^^^£,2^-3)-2/75(2-2a)/7[Qj^ ^ 3 _ 

where Di is the luminosity distance in Mpc, a is the spectral index, 8 is the Doppler factor, Vm 
is the turnover frequency for synchrotron emission in GHz, Sm is the flux density at the turnover 
frequency in Jy and R is the linear radius of the emitting region in cm at the turnover frequency. 

3.3 Distance to SMBH 

The relativistic jet model was first proposed to explain the variability of non-thermal emission 
in quasars and blazars [5, 6]. As Marscher et al. [14] puts it: 

The jet is assumed to be generated at some point /?apex. beyond which it flows at a 
constant Lorentz factor L (at speed jSc), confined to a cone of constant opening half¬ 
angle 

If the jet is confined only by its own inertia, the jet density will decrease as 1/ r^, where r is the radial 
distance from the jet apex. If the magnetic field is parallel to the jet (toroidal), the magnetic field 
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will decrease as B|| oc when it is perpendicular (poloidal) oc This can be re-expressed 
as B oc ( 1 /r)” [e.g. 10 ], where n is an exponent for toroidal (n=l) or poloidal (n= 2 ) magnetic field 
configurations. We can compute an estimate of the distance from the “core” to the jet apex with the 
ratio of magnetic fields and the separation between the “core” (Cl) and downstream component 
(C2): 


Bc\ 

Bc2 


rc\ +Arc2 
rc\ 


(3.2) 


Where rci is the distance to the black hole from the “core” in mas, Arc 2 is the separation between 
the “core” and stationary feature. Rearranging, we get: 


rci 


Arc2 

{Bc\/BciY — 1 


(3.3) 


This can be generalised to be able to solve for the magnetic field at any arbitrary location along the 
jet: 


B2 _ TBil 

Bi R2 


(3.4) 


where Bi and B 2 are the magnetic fields af arbitrary transverse radii Bi and R 2 . Hence if the “core” 
is assumed to be resolved, we can estimate the magnetic field sfrength at a distance at any arbitrary 
distance (e.g. IOB 5 ). 


4. Discussion and Conclusions 


1.4 
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Figure 1: Equipartition magnetic fields, binned by de-projected distance from the “core”. Red bars indicate 
the average upper limits and yellow bars indicate the average lower limit. The error bars indicate the scatter 
of the limits. 

Although only upper and lower limits on the magnetic field strength can be computed, we 
can bin and average these limits. These averaged magnetic field limits are plotted in Fig. 1 as a 
function of de-projected linear “core” separation and over all sources. De-projected distances are 
binned ever 2 parsecs over the first 10 pc. From this, we find fhat the Doppler corrected magnetic 
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field strength in the “core” is on average between ~0.8 and 1.4 G, which is stronger than the re¬ 
sults found via “core-shift” by O’Sullivan & Gabuzda [19], which on a sample of 6 quasars found 
magnetic field strengths of 10s to 100s of mG. These values are not necessarily inconsistent with 
our results as these values were measured at cm-wavelengths and hence further downstream than 
the mm-wave “core”. Our values are however consistent with the values found in the MOJAVE 
sample by Pushkarev et al. [20], which found values of ~0.9 G in FSRQs and ~ 0.4 G in BE Eacs. 
In our small sample of six sources, we find local maxima in the B-tield distribution which may 
be interpreted as possibly due to recollimation and oblique shocks [e.g. 13], but due to the low 
sample size, we cannot rule out any scenario. We can also test the rate at which the magnetic field 
decreases with increasing “core” separation, by performing a power-law fit to the limits, B oc 
This yields a value ofn = 0.3 ± 0.2. 

Using the decrease in magnetic field strength of 0.9< B <1.4 at zero “core” separation to 
0.2< B <0.7 at a de-projected separation of 4 —6 pc, yields an estimated distance from the jet apex 
to the mm-wave “core” of between ~ 1-3 pc. This likely places the mm-wave “core” outside of 
the broad line region. Taking a typical SMBH mass of Mbh = 5x 10^ solar masses and assuming 
a toroidal (n=l) geometry, we can extrapolate the average jet magnetic field strengths back to a 
distance of IOB 5 , yielding strengths of Bapex ~ 5 — 20kG. This points towards a magnetic launch¬ 
ing mechanism of the VEBI jets, as recently described in the MAD and other similar models [e.g. 
16, 17, 24] and is near the maximum Eddington magnetic field [e.g. 21]. An observational study 
by Silant’ev et al. [23] also found magnetic fields of the order thousands of Gauss in AGN, inde¬ 
pendently confirming such high values. 

Zamaninasab et al. [25] showed that accretion disk luminosities are related to their black hole 
masses and hence to their magnetic field strengths. Eor black-hole masses of > 10^ solar masses 
and higher accretion disk luminosities, the magnetic field estimates can be Biorj ~ 30 kG or higher, 
consistent with the estimates derived here. However, if the magnetic field configuration is consis¬ 
tent with n> magnetic fields would be in the MG range, considerably stronger than expected 
from MAD and higher than theoretical Eddington limits, supporting the view that blazars have 
toroidal magnetic field configurations. 

5. Conclusions 

We have presented a novel approach to determining the magnetic field geometry of AGN jets. 
By spectrally decomposing 6 AGN using 3 mm and 7 mm VEBI, we have determined limits on the 
average magnetic field strength gradient. We have found that the average magnetic field strength in 
the mm-wave “core” is between ~0.8 and 1.4 G, decreasing to 0.3-0.4G at a de-projected “core” 
separation of 3-5 pc. Fitting a power-law to these values (B oc r^”), yields n = 0.3 ± 0.2. By 
assuming a conically expanding jet and a toroidal (n=l) magnetic field geometry, we estimate the 
jet apex to be ~ 1-3 pc upstream of the mm-wave “core” and likely outside of the broad line region. 
Extrapolating the limits on the magnetic field strength to near the apex of the jet suggests magnetic 
field strengths of Bapex ~ 5 — 20 kG. We find that such strong magnetic fields are broadly consistent 
with magnetic jet launching, for example under the scenario of magnetically arrested discs. 
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